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Abstract 

The rate of emission of heavy dileptons from QGP is found without an assumption 
of its complete thermal equilibrium. We base on the real-time quantum field kinetic 
approach [1] and use the expansion up to the second order with respect to strong 
coupling constant g. The final answer is not free from the collinear singularities and 
we show that this is the actual issue. As a result the main contribution to the rate 
of the heavy dileptons production at M/T ~ 10 comes from the process qqg — > 7*. 
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Introduction 

Transparency of hadronic plasma for electromagnetic signals inspires hope that 
they may carry information about the intimate details of interactions in the quark- 
gluon media. For this reason they are considered as promising probes of the hot 
hadronic matter [2]. Initial theoretical efforts were put basically to the study of an 
influence of hydrodynamic background on the total yield of the dileptons [3]. The 
production rate was taken from the first Born's diagram, qq — > 7*. The dynamics 
of quark-gluon interactions was present only virtually via an assumption that it has 
led to a thermal distribution of quarks. The further improvement of the microscopic 
picture also took into account the gluon interactions [4] up to the a s -order. 

The most encouraging result was obtained by Baier, Pire an Schiff [5]. They 
found that the rate of emission of the dileptons in the a s order is defined as well as 
the first Born's term: all the infrared and mass singularities has cancelled out in their 
approach. It was the direct check of the Lee-Nauenberg theorem for the emission of 
dileptons from the equilibrium thermal bath. This result was confirmed in somewhat 
different manner in [6]. 

The most important consequence of this study was a statement that even apart 
from the factor exp(— M/T) the rate of dilepton emission is a decreasing function of 
M/T. 

Nevertheless there remain some questions. First, any deviation from the detail 
balance relationships of the true thermal ensemble breaks a zero balance of mass 
singularities. The calculation given below is at least an example. Second, the rate 
of dilepton emission in the process qqg — > 7* increases (apart from exp(— M/T) ) 
with the growth of the dilepton mass. The effect is due to specific features of the 
phase-space volume of the reaction and is not sensitive to the nature of the infrared 
cut-off. The leading role of this reaction was noticed in [7] . Recently it was confirmed 
[8] by numerical calculations. In this paper we give an analytic estimate. 

We confine ourselves to the same first a s -order also. The main difference of this 
work from the previous ones is that we consider the quark gluon plasma without 
thermal and chemical equilibrium. We show that both effects, the absence of collinear 
singularities and suppression of the process qqg — > 7*, are the artifacts of ideal thermal 
equilibrium. If the equilibrium is broken then at M/T ~ 10 the rate of this process 
overwhelms the Born's rate while the rates of all the other processes are suppressed 
at least by the factor T/M. 

An initial goal of this study was to design a tool which could have been used for 
calculations of the dilepton emission against the background of partons distributions 
generated from cascade simulation of the A — A collision. Some scenario show [9] that 
cascade do not reach thermal equilibrium and that quarks are relatively suppressed 
with respect to hot glue. So we used some flexible analytic approximation for dis- 
tribution of the partons in favour of the opportunity to get an analytic answer and 
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reliable estimate of the leading terms. 

While choosing this approximation we kept in mind the following scenario of the 
heavy ions collision. The initial stage of a collision at RHIC or LHC energies (r ~ 
0.5 fm) is a region of nucleons fragmentation and development of the initial parton 
cascade. Dileptons are emitted only due to the hard Drell-Yan process. Our region 
begins a little bit later when partons are already chaotized and may be described by 
the one-particle distributions. 

The most general density matrix which simulate any given ahead form of the 
one-particle distribution is of the next form, 

p = e ~fj {N,P)a+ (N,p) aj (N,p) ^ j ^ 

N p,j 

where N label the space cells on the hypersurface of the initial data and rij(N,p) = 
a~j~(N,p)aj(N,p) is an operator of the number of partons of the sort j and quantum 
numbers p in the N-th sell. Thus we completely neglect all the correlation effects in 
the phase space of the partons. 

The density matrix (1) allows one to apply the field theory to the description 
of the further evolution and results in a set of habitual elements like cross-sections, 
self-energies, vertices. If the functions fj(N,p) which define the shapes of partons 
distributions are completely arbitrary then we are totally confined to numerical cal- 
culations. The only piece of theoretical study will be connected with the necessity to 
transform some singular expressions to the form when numerical calculations will be 
unambiguous. 

Allover this paper we use the Boltzmann-like distributions damaged by intro- 
duction of specific parameters (q and (g which all together will be considered as the 
measure of nonequilibrium in the quark-gluon system. We adopt for quarks and gluons 

n F (p) = CQe HH/T , n B (p) = CceHH/T (2) 

where u**, the 4-velocity of the nonequilibrium partonic media, fugacities ( and tem- 
perature are very smooth functions of space-time coordinates. Scenario of the hot glue 
[9] leads to (q < (g < 1 an d T is rather a formal parameter than a thermodynamic 
temperature. We hope that distributions (2) are enough flexible to approximate the 
cascade approaching to the thermal equilibration. 

An important advantage of this parametrization is that all the calculations can be 
performed analytically to the very end and result in no more than one-dimensional 
convergent integrals. 

The paper is organized as follows. In Sec.l we briefly remind the basic definitions 
as they were derived in ref. [1] and trace the way of perturbative expansion assuming 
the weak coupling between the quarks and gluons. In Sec. 2 we split the general 
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expression for the rate of the dilepton emission into the pieces which at some heuristic 
level may be called as real processes and radiative corrections. In Sec. 3 the rate of 
the real processes is divided into more specified parts, annihilation and Compton and 
the integration over the distributions of initial partons is carried out. We carefully 
analyze the kinematic regions where the further radiative corrections will be needed. 
In Sec. 4 we do the same for the virtual corrections. In Sec. 5 we perform the remaining 
integration over the final state of the partons. At this step we still have a set of the IR- 
divergent one-dimensional integrals. Sec. 6 describes the assembling procedure which 
results in the IR-finite answer for the rate of emission. We show that scale which is 
responsible for IR-finiteness of the observable rate is T (not gT\) and put forward 
the arguments why the double logarithms should be exponentiated into the K-factor. 
In conclusion we analyze the mass singularities and show that they are naturally 
restricted by the amplitude of the forward scattering of hard quark in the partonic 
media. By some chance this cut-off coincide with the so-called thermal mass. This 
mass can not be associated with the rest mass of any quasi-particle. 

l.Perturbative expansion for the rate of emission 

For sake of completeness we shall start with reminding the definitions and the main 
framework of calculations [1]. We shall close this section by the formal expansion of 
the rate of emission defined in the scheme of quantum field kinetics up to the first 
order in a s . 

The inclusive rate of the dilepton emission is given by 



where k = ki + k 2 , L^ v = k!{k v 2 + k^k\ — g^ u (k 1 k 2 — ml), is a trace of lepton spinors, 



is a kind of photon Wightman function averaged with pQ CD and e is an electric 
charge of the lepton. 

We suppose that an explicit separation of long-range and short-range orders takes 
place and have passed from the inclusive cross-sections to the emission rates per unit 
volume. That is we have assumed that even not reaching thermal equilibrium cascade 
produces at least locally homogeneous distributions of the partons. 

The integral equation for A io can be effectively solved in a manner described 
in [1] and we get the rate of dileptons emission expressed via polarization operator 




(1. 1) 



Aft(-fc) = -ij d\x-y)(T + {A^x)S + )T{A»{y)S))e 



ik(x-y) 



(1. 2) 



nio(-fc), 
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,o,o dN e+e - _ . 2 L IIU (k 1 ,k 2 ) Tl^(-k) 



where [k 2 ]~ 2 stands for a product D ret (k)D a d v (k) out of the photon mass shell. 

Integrating this distribution over the relative momentum of the leptons we get a 
distribution of the dileptons over their masses and total momenta: 

dN e+e - = dN e+e - iirelg^U^j-k) 

d A kd 4 x °d 3 kdM 2 d 4 x 3(2tt) 6 M 2 1 ' } 

Polarization matrix in these equations is given by an expression 
IWs,i/) = -t(-l) A+B e 2 £ (-l) R + s [dZd V YG AR (x,OE» mB (Z, m y)G S A(v,x), 

R,S=0 J 

(1. 5) 

where e stands for the electric charge of a quark and a subsequent summation over 
the u- and rf-quarks is assumed, e 2 = (5/18)e(j. It contains an electromagnetic vertex 
dressed by strong interaction, 

^(».v ; ») = (-i)^ +p f[G :' (l : ,) . lm . (i.6) 

The matrix Schwinger- Dyson equation can be solved in the form, 

[G-'Ub = [Go V - M AB , (1. 7) 

where the matrix of quark self-energy looks as 

M AB (x,y) = i{-l) A+B gl £ {-l) R+s I dZd V x (1. 8) 

R,S=0 J 

and T is the vertex of strong interaction. The inverse Greenian Gq 1 is nothing but 
the Dirac operator with the "external" field A{x), 

[G^\ AB {x, y) = (-l) A S AB S(x - y)[id x + etA(x A ) - m] (1. 9) 

Proceeding from the Eq.(1.6 ) and confining ourselves to the two lowest orders of the 
vertex expansion we get 



nft(-fc) = te 2 N c J -^l"G 10 (p - k)YG 01 (p) + e 2 g 2 N c C F EJ" 1 )* 4 " 5 x 



i?,S=0 



x /^^ J D5 R (g)7 M G 1R (p-A;)7 A G R o(p + g-A;)7 1 'Go5b + g)7AG sl (p)(l. 10) 



The second term of this equation is already restricted to the a s order. The further 
expansion of the exact Greenians in the first term is performed by means of equa- 
tions^], 

Goi,io = Coi,io — G ret M i t i G a dv + G rei M ret Goi,io + Goi t ioM a d v G a( iv, (1. 11) 

where in the order we are interested in all the propagators are to be taken bare and 
self-energies should be calculated to the first nonvanishing order in strong coupling 
constant. The later are as following, 

M AB (s) = *(-l) A+B |^ / dSl,G AB {s + qh x D x / A (q), (1. 12) 

and 

M retMv (s) = J d A q[-f^G retMv (s + qhxD^(q) + ^G^s + qh\D^ v ret (q)}. 

(1. 13) 

2. Separation of real processes and radiative corrections. 

Eqs. (1.11) - (1.14) give us the rate of emission in the closed form containing 
all contributions of all possible processes. To split them we should notice that an 
indication of the type of the process comes from its statistical weight. Thus we 
can single out the direct annihilation, qq — > 7*, annihilation with the emission or 
adsorbtion of a gluon, qq — > #7* and qqg — > 7*, and Compton scattering of quarks on 
a gluon, qg — > qj* and qg — > qj* . 

If in the first line of Eq.(l.ll) we substitute instead of exact Greenians Gio,oi their 
unperturbed values we get the first Born's term, related to direct process qq — > 7*: 

^Bom(-k) = ie 2 N c j -0^rG lo ( P - k)YG 01 (p). (2. 1) 

The diagram of the Born's term is given at Figl. 

If we insert terms with the off-diagonal elements of the self-energy we get two 
equal terms related to the real processes, which contribute to the second Born's 
approximation term 

Iin-*) = 2eYN c C F J^-D w (q)Tr{rG w (p-k)Yx 

xG ret (ph x G 01 (p + qhxG adv (p)} (2. 2) 

The second contribution to the real processes comes from the two equal items with 
R 7^ S in the sum in Eq.(l.ll). 
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ICC-*) = -2e*g*N c C F f d j0D w ( q )x 

xTr{^G w {p - kh x G 00 (p + q - k)>fG 01 (p + q)lxG u (p)} (2. 3) 

In a short while we shall rearrange H re ai — n a +IIb into a sum Tl rea i = n em +n a b s +Il com 
of the inclusive subprocesses of the gg-annihilation with the emission or absorption of 
the gluons and the Compton (bremsstrahlung) process. The diagrams of II a and U b 
are given at Fig. 2. The crossed lines present the on-mass-shell partons with density 
matrix (1). Their statistical weights originate from the off-diagonal densities Gab 
and Dab, A B, which carry a certain sign of the energy. 

The rest items from equations (1.11) and (1.12) are the self-energy and vertex 
corrections. 

The two remaining terms of the Eq.(1.12), containing M ret and M adv give rise 
to the radiative corrections of the self-energy type. They account both for the vac- 
uum effects and for the amplitude of the forward scattering of quark on the partons 
produced by the cascade. 



^aJ-k) = e 2 g 2 N c C F J ^0Tr{^G w {p - k)-fG i(p) x 

x{[l X G ret (p + q^xD^q) + ^dip + q) lx D adv (q)}G ret (p) + 

+ h X G adv (p + qhxD 1 (q) + ^G^p + qhxD ret {q)}G adv (p)} (2. 4) 

This expression contains several functions with overlapping singularities and it seems 
better to transform it 

Wnassi-k) = eV^||0Tr{ 7 A G s (p)7^ lo (p-A;) 7 ^o 1 (p)7A}x 

x(G s (p + q)D 1 (q) + G 1 (p + q)D s (q)) (2- 5) 

The items with R = S in Eq.(l.ll) produce the vertex-type radiative corrections, 

d 4 pd A q 



KU-k) = - e2 9 2 NcC F J^0x 



(2nY 

xTr{YG 10 (p - A;)7 A G o(p + q - k)YG 00 (p + qhxG 01 (p)D 00 (q) + 
+YG w (p-kh x G 11 (p + q-k)YG n (p + q) 1 xG i(p)D n (q)} (2. 6) 

The vacuum part of the 7r ver t is reasonable to calculate just in the form given above. 
It is due to the T-independent parts of Goo and Gn- For the thermal part it is better 
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to transform, 

KU-k) = -e 2 g 2 N c C F J ^|Tr{ 7 A G 01 (p)7^ 10 (p - k)^ x 

x (G s (p + q - k) lv G s (p + q)D 1 (q) +2G 1 (p + q- h)>fG a (p + q)D 9 (q)) (2. 7) 

These transformation present "real" and "virtual" processes in the same form: now 
all thermal diagrams have three cuts! The difference is that in radiative corrections 
the extra cut correspond to "elastic scattering on the distribution of the partons ". 
The typical graphs of this corrections are plotted at Fig.3. Statistical weights of the 
additional cuts always originate from the local density of states G\ (or D^.They are 
contributed both by induced emission and absorption. 

3. Real processes. Integration over the initial states. 

In this section we begin calculation of the real processes. We do not assume that 
quarks and gluons are in thermal equilibrium and use the partons distributions (2) 
for the analytic calculations. 

At the end of this section the reader will also find physical discussion, what part 
of the phase-space is important for the heavy dileptons production. 

In what follows we shall calculate only the distribution of the emitted dileptons 
over their total 4-momenta, so we calculate only n(k) = g^TV^iy—k). 

For the Born's term the first integration is the last one, 



/,x 3ie 2 . 2m 2 . I 4m 2 ,,, , /T 
n Born (k) = -^(1 + ^1 ~ ¥ MV^. (3. 1) 

In the next perturbative order we have, 

Kreal = ~ ^^ L J d^q5(q 2 )5[(p - kf - m 2 ]5[(p + qf] X 

x (SW em + SW abs + SW com ){l + - A + - - - + ° } (3. 2) 

[p z — m z y p z — m z p z — m z + Ikq 

where we denoted: 

A = 2m 2 (k 2 + 2m 2 ), B = 3m 2 + k 2 + 2(kq) - — ~ ^ 

2{kq) 

The expression in the curly brackets is nothing but a sum of the squared moduli 
of the matrix elements of the annihilation processes with emission, qq — > g^y*, and ad- 
sorbtion, qqg — > 7*, of a real gluon or the Compton process, qg — > qj* and qg — > qj*. 
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Just for the purpose of the following calculations it is written not it terms of habit- 
ual Mandelstam variables (s,t,u). Specification of the process is due to statistical 
weights. 

The statistical weight of the annihilation process with emission of a gluon looks 

as 

SW em = 9(k - p )9(q + p )9(q )n F (k - p )n F (q + p )[l + n B (q )] « 

» 9(k - p )9(q + p )9(qo)C 2 Q e- {ku)/T (e-^ )/T + Cce" 2 ^) (3. 4) 

The statistical weight of the annihilation process with absorption of a gluon is 

SW abs = 9(k - p )9(q + Po)d(-q )n F (k - p )n F (q + p )n B (-q ) ~ 

» 0(k o - p )9(q + Po)0(-qoK 2 Q CGe- {ku)/T . (3. 5) 

For the Compton rate of the dilepton emission the statistical weight equals to 

SW com = -0(p o + qo)9(p - k o )0(-q o )n F (p o + q )[l - n F (p - k )]n B (-qo) - 
-0(-po - qo)9(k - p )9(-q )n F (k - p )[l - n F (-p - q )]n B (-q )- « 
« -CqCg{0( P o + qo)0( Po ~ k )9(-q )e-^ T [l - CQ e-^- ku V T } + 
+0(-Po ~ qo)0(ko - Po)9{-q Q )e' ku / T e^ u+ ^ T [l - C Q e^ + ^/ T ]}(3. 6) 

The 4-vector k^ + q^ is time-like. So we can perform integration over p using 
the Breit reference system where k + q = 0. The exact analytic result has the next 
covariant form: 

vw = ^^f^ le-^ j rf 4 g5(g 2 )^(g )^(A: + g) 2 -4m 2 ]^(A;g){e^ + C G e- 2 ^ T }, 

(3. 7) 

Tata = ff J 4 C F ( 2 Q ( G e- ku/T J d*q5(q 2 )9(q )9[(k - qf - 4m^ a (-kq) (3. 8) 
where 



M 2 + 2m 2 M 4 -4m 4 1 + V 1 " 4m 2 /(M 2 + 2x) 

Fa{x) = -(1 + + — = ) In v + 

^ 2a; 1 - ^/l - 4m 2 / (M 2 + 2x) 



,„ M 2 + 2m 2 M 2 (M 2 + 2m 2 ), / 4m 2 , n ft , 

We emphasize that the " nonequilibrium" parameters £q and Cg carry a significant 
additional information. They allow to separate contributions of spontaneous and 
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induced processes. This information would have been important even in the case 
of true thermal equilibrium but it is completely hidden if we use the detail balance 
relations at the early stage of calculations. For instance, examination of the expression 
for ir em immediately shows that the first of the "gluon exponents" originates from 
quark and anti-quark distributions and momentum conservation. It describes the 
spontaneous emission of a gluon and it is not sensitive to a shape of the initial gluons 
distribution. 

Practically, introduction of £'s prevents the eventual cancellation of mass singu- 
larities between real processes and vacuum or thermal mass and vertex radiative 
corrections. The level of vacuum fluctuations and the rate of spontaneous processes 
are fixed de fault by ( vac = 1. 

For the Compton rate we start with a chain of changes of variables aimed to reduce 
the statistical weights to the uniform shape. These are p — > — p + k — q, in the first 
term and q — > — q — p in both terms. Then we may perform an exact integration over 
p using the same Breit reference system. It gives 

W = - ^^f' CQCGe-^ j d*q5(q 2 - m 2 )6(q )F c (kq) [e^ - CqC 2 ^] 

(3- 10) 

where q is the momentum of the (anti) quark in the final state and 



. . 4m 2 + M 2 - 2x - 2(M 4 - 4m 4 )/(M 2 + 2x) , m 2 + x + ^x 2 - M 2 m 2 

J- C {x) = In 

2\/x 2 — M 2 m 2 m? + x — V x 2 — M 2 m 2 

4(M 2 + 2m 2 ) (M 2 + 2x){m 2 + x) 

+ M 2 + 2x + (M 2 + m 2 + 2x) 2 ^ ' ' 

Besides an opportunity of analytic integration another advantage of the Breit system 
is that it reveals some additional information about the virtual state of the quark 
field. Namely, in annihilation process the energy p of virtual state in the Breit 
system equals to ratio of squared dilepton mass to twice energy of initial quarks (or 
dilepton plus gluon). For the most profitable configuration of heavy dilepton and 
soft gluon this is about half of dilepton mass; the 4-vector p^ is space-like and lies 
very close to the light cone. For the same reasons the energies of the initial quarks 
are almost the same but their 4-momenta are slightly time-like. For the Compton 
process we get the same estimate. In the most probable case of heavy dilepton with 
low momentum k the Breit system almost coincide with the rest frame of the media. 
So when we will need further mass corrections to the quark propagators they should 
be calculated just at these momenta, p ~ M/2 and p 2 ~ 0. 

Another important information is that at any finite variable q the very integrands 
in Eqs. (3.7,8,10) are divergent at m — > 0. This mean that the collinear divergen- 
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cies have explicitly different origin than the IR-divergencies and should be treated 
separately. 

4. Radiative corrections. Integration over the initial states. 

This section is technical. It could have been safely berried under the common "As 
it can be shown..." if there were no difficulties associated with the extremely singular 
form of the integrand. Being worked out naively they can essentially change the final 
answer. 

Following Eqs.(2.5) and (2.7) we calculate contributions of quarks and gluons from 
the initial partons distributions to the vertex and self-energy corrections separately. 

The vertex radiative correction is divided into the two parts contributed by the 
Fermi and Bose parts of the partons distribution, U vert = n^ ert + U.^ ert : 

< rt = %e2 ^f F J d*pd*q5(q 2 )5[(p - kf - m 2 }6[(p 2 - m 2 )]SW v B ert X 

4m 4 - k 4 - 2kq(m 2 + k 2 ) - 4(kp) 2 + 4pq(k 2 + m 2 + kq) 
X ((p + q - k) 2 - m 2 ){{p + q) 2 - m 2 ) ' ( ' ' 

<ert = %(?9 ^f F I d*pd 4 q6(q 2 - m 2 )6[(p - k) 2 - m 2 ]6[(p 2 - m 2 )]SW v F ert X 



-m 2 k 2 - 2kqm 2 - 4(kp) 2 + 2pq(k 2 + Am 2 + 2kq) 
(q + k) 2 - m 2 )(p - q - k) 2 



In the same way we single out the two parts of the mass corrections, IT r 



IF + IF : 

mass 1 mass 



b __ ie 2 g 2 N c C F r 4 



^mass 



J d 4 pd 4 q6(q 2 )5[(p - k) 2 - m 2 ]5[(p 2 - m 2 )]SW* ass x 



V V 4 

(p 2 — m 2 ) ((p + q) 2 — m 2 ) 



<ass = ^^f' J d 4 pd 4 q5(q 2 )5[(p - k) 2 - m 2 ]5[(p 2 - m 2 )]SWf 



V 

mass 



V V 

(p 2 — m 2 ) (p — q) 2 ^ ^ 

where V means the principal value of the integrals and T1Z stands for the trace of 
spinors from the corresponding Greenians, 

T1Z = 2[{p 2 — m 2 ) 2 + {kq + pq — kp — m 2 ){p 2 — m 2 ) — 2m 2 (pq — kp) — 2{kp){pq) + 2m 4 ] 

(4-5) 
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The corresponding statistical weights for these graphs are 

SWvert,mass = 6 ( k o ~ Po)0(Po)n F (k - p )n F (p ) [1 + 2n B (\q \)} « 

« 9(k - p )9(p )C 2 Q e'^/ T [l + 2( G e-^ T ] (4. 6) 



SW£rt,mass = 0(k - p )9(p )n F (k - p )n F (p )[l - 2n F (\q \)} w 

« ^(A; -Po)e(po)( 2 Q e~^ T [l - 2^-1^]. (4. 7) 

One may notice that two of three statistical factors are modified by the theta- 
functions. These are due to the real partons in the initial or final states. The third 
factor works at both signs of energy which means a simultaneous account for both 
absorption of quanta from some state in the partonic bath and backward emission to 
the same state. 

The first integration d A p over the initial partons distributions is naturally per- 
formed in the reference frame of the dilepton. For the vertex corrections they result 
in 



b _ *e 2 g 2 N c C F 

n vert 47f 4 



( 2 Q 9(k 2 - Am 2 )e- ku / T J d*q5(q 2 ){\ + C G e~^ T } x 



i 4 -4m 4 , 1+v . x 



< ert = lC 9 *f F e Q 0{k 2 - 4m 2 )e-^ T x 
x / A% 2 - m 2 ){ 1 - - C Q e-^ T }[ k \ 2 k l~ k 2 q m2 v + (4. 9) 



| (k 2 + 2m 2 )(kq + m 2 ) - m 2 (k 2 + 2kq)/2 ^ 2m 2 + kq + v^](kq) 2 + k 2 m 2 



(k 2 + 2kq)y (kq) 2 + k 2 m 2 2m 2 + kq — (kq) 2 + k 2 m 2 

When calculating the self-energy corrections we meet several ill defined products 
of the singular functions like x5(x)V(l/x). They need much care and the way to 
handle them is described in the Appendix A. This is hardly the most cumbersome 
part of calculations. It results in 



7T 



Lss = ^^f^ ^k 2 - Am 2 )e^J d\5(q 2 ){\ + ( G e^} x 



m 2 (fc 2 -2m 2 ) 1 , l + v ,2m 2 k 2 + 2m 2 k 2 (k 2 + 2m 2 ) 1 _ . A 

x « 1+ Wd "°— + — k, W^ w (4 ' 10) 
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rLss = ^^f' Cle^ - 4m 2 )e"^ J d * q5(q 2 -m 2 ){\- C Q e~^} 



r k 2 + m 2 + kq , 2m 2 + /eg + vJikq^Tk 2 ^ 2 2m 2 . . 

x[ =ln / = ~ T2~ u ]» ( 4 - U ) 

J{kq) 2 + k 2 m 2 2m 2 + kq — vJ (kq) 2 + k 2 m 2 fc 

where we denoted, v = \Jl — Am 2 /k 2 . 

Terms corresponding to 1/2 in curly brackets originate from vacuum fluctuations. 
They are both IR- and UV-divergent. Removal of this divergence is a subject for 
renormalization. We adopt a standard procedure of the on-mass-shell renormalization 
in the asymptotic states both for quark self-energy and electromagnetic vertex [10]. 
Thus we insist that vacuum part of the quark self-energy M(p) has a zero of the 
second order at p 2 = m 2 . In the same way we adopt that when all three momenta of 
the vertex function, 

n^k 

T»(k 2 ) = ^F^k 2 ) + —^F 2 (k 2 ), (4. 12) 

are on the asymptotic mass shells of the free particles then the form- factor iq(0) = 0. 

This kind of renormalization immediately wash out vacuum part of the quark self- 
energy and leaves the following contribution of quark electromagnetic form-factors to 
the dilepton emission as the residue [10]: 

7r™ c t = ie2a 2f F C Q 0(k 2 - Am 2 )e- ku ' T M 2 x (4. 13) 

x{-Kv) [ ^+ (1 + %> - [f + £v] In g + 2(1 - jr ihv ^} 

The IR-divergent integral dq/q originates from the soft space-like gluon exchange 
between the fermion legs of the electromagnetic vertex and the function h(v) is given 
by 

i / \ 4m 4 . 1 + v 2m 2 . , A . 

%) = -(l-^)ln T ^ + (l + ^ (4. 14) 

An asymptotic of the last integral in (4.14) at m/M « 1 looks as 



/' 

Jo 



atanhv xdx ^ \ \ _|_ v \ M 2 tl 

w h - In 2 S3 - In 2 — (4. 15) 

tanh xG4 1- v4m 2 



and produce large negative contribution to the total rate of emission. This double 
logarithm appeared from the collinear configuration of the vertex and comes to be 
infinite for massless quarks. It may be quite safe if it is properly balanced with 
the terms of the opposite sign produced by the other processes which are coherent 
with this one. As it was shown in [5,6] such a tiny balance takes place in the ideal 
thermal ensemble. For our distributions it does not happen. But after examination 
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of the scales of different processes we shall see that this double logarithm should be 
exponentiated into the K-factor. 

5. Final integration. 

There are two time-like 4-vectors in the theory, the 4-momentum of the dilepton, 
k^, and the 4- velocity of the media, u^. So we may continue calculations either in the 
rest frame of the media or in the rest frame of the dilepton. We will confine ourselves 
to the rest frame of the media, u = 0, basically because it is more natural to work 
with the isotropic distribution. 

Let us consider the term responsible for the gg-annihilation accompanied by spon- 
taneous emission of the gluon. Its "weight" is CqCvoc = Cq- After using of the delta- 
function and one trivial asimutal integration we come to the two-dimensional integral, 

-ttT J d[e^ T ] J^dzqF a (k q-\k\qz). (5. 1) 

Here and further on k stands for ku and |k| 2 = (ku) 2 — M 2 . The trick which allows 
one to reduce this double integral to a simple quadrature is as follows [11,12]: 

i) change the angular variable z = cos 6? for the new one, q(ko — \\t\z) = y; 

ii) integration with respect to q by parts: in a miraculous way the integral dy 
disappears ; 

iii) changes of variables: q = k + y (or q = k_y) where k± = k ± |k|. 
The result reads as 

— {lim / (1 - e-y/ k - T )F a {y)dy + / (e^ k - T - e-y^ T )T a (y)dy} (5. 2) 



hm r + ^(/(0)+yf(0) + ...) = /(0) In |++0(A), (5. 3) 



Though both lower and upper limits of the first integral tend to zero when A — > the 
very integral is finite. Indeed, if f(0) ^ then 

rXk+ dy , i fl( ^ i N x/n ^ k^ 

y 

and 

r Xk + dy , 2|k| 

III A I — ' 1 ' ' " 

A^o J\ k _ y 1 

Foreseeing the future cancellation of soft divergencies we can safely put the lower 
limit equal to zero and wright the final result for the annihilation with the gluon 
emission in the following form, 

roo poo 

+2/ S{y)F a {M 2 y)ydy + 2( G / S{2y)F a {M 2 y)ydy\. (5.5) 

JO JO 
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The terms with the extra factor (g relate to the induced emission of a gluon and we 
denoted, 

s w = e ^W (5 ' 6) 

In a quite similar way we obtain the expression for annihilation with the absorption 
of a gluon, 

= *^£=<a&e-TM«[(-l + ^ h |± )M .„) + 

r(l-4m 2 /M 2 )/2 



+2 / F a {-M 2 y)ydy\ (5. 7) 



where it is helpful to denote v(y) — y 1 — 4m 2 / (M 2 + 2y), (now t>(0) = i>) and rewrite, 
2 N 2m 2 1 { 4m 4 1 l + v(y) 



T a {M 2 y) = -(1 + (1 + — ) - + (1 - — T )— ) In 



M 2 'y v M AJ 2y 2J 1 - v(y) 
2m 2 , 1 . 2m 2 , 1 , . , 

+(1 + < 1+ m^ + (1+ aF>v W! '> (5 ' 8) 

Assembling all gluon contributions to self-energy and vertex radiative corrections 
together we get 



POO 

+2 / (K(y)+K(-y))S(y)ydy] (5. 9) 

Jo 

where function lZ(y) in the integrand has a form 

Equations (5.5), (5.7), (2.9) and the vacuum contribution to the vertex (4.13) form 
a closed system. Though each of them contains IR-divergent integrals over gluon 
momentum, when taken together, they produce finite result. 

Further calculation of Compton process and radiative corrections from fermion 
part of the initial state background meet no infrared problems at low momenta of 
gluons. So we can safely choose the simplest way of calculations using the rest frame 
of a dilepton. 

7W = - e2g l^ CF CQCGe- ku/T 2m'f -0=F c (x)[C(x,T) - C Q C(x, |)] 

(5. 11) 
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where 

Mlklm \ 
MT_ X ) 

-x) 

MT ' 



c(«.n = «*-5SSn^ ,5.12) 



and 



_ 4(M 2 + 2m 2 ) (M 2 + 2Mma;o)(m 2 + Mmx ) 

cW ~ M 2 + 2Mmx + (M 2 + m 2 + 2Mmx ) 2 ( j 

1 , „ „ , 2 rt „ , 2(M 4 - 4m 4 ) . , m 2 + Mm(a;o + x) 

-{4m + M — 2Mmx — — ^ — — — -} In 



2Mmx L M 2 + 2Mmx J m 2 + Mm(x - x) 

where xq = \/l + x 2 . The contribution of the Fermi part of intial partons to radiative 
corrections is of the next form: 



^ rad 



- lC 9 ^f F ^ T 2mM jT JL*L.(j: f ( x , Xq ) - f f ( x , -x ))C(x, T) 

(5. 14) 

where 

^ ( "' Xo) = ^ ( M 2 + 2Mmx 2H5 ' 15) 

+ 1 ( (M2 |'r 2)( ^ + mMX ° ) - ^ - - 2 - -^o) In ^ + Mm(*o + *) 
x v (M 2 + 2Mmx ) 2 y 2m 2 + Mm(x -i) 



6. Assembling and cancellation of the IR-singularities. 

The final answer for the rate of emission is given now as a sum of many terms, 



dN e+e - me 2 

:7l(k), (6. 1) 



d 4 kd*x 3(2tt) 6 M 2 

n(k) — IT Born + '"em + '"abs + '"'com + '"fad + 71 'rad + 71 ~ 



vac 
verti 



which are given by Eqs. (3.2), (5.7), (5.10), (5.15) and (4.13). Most of them are the 
IR-divergent one-dimensional integrals. Nevertheless the whole sum is IR-finite. We 
show now how it happens and find out the physical scale which is responsible for this 
phenomenon. 

The cancellation of the infrared divergencies due to the soft gluons should be 
looked at separately for terms with the different powers of fugacities ( G and (q. The 
greatest contribution to the rate of the emission of heavy dileptons comes from the 
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annihilation with the absorption of a gluon. The sum of the leading IR-singular terms 
from 7r a b s and n ra d results in the finite integral, 

n abs « ^(^( G e- h /V f {-(1 - §r)(L(-<z) - e^L(O)) + 

+(l + ^)(^(-?)-e-^(0))} (6. 2) 

where q max = (M 2 - 4m 2 ) /2M ~ M/2, and 



•w-^-iPTW ^-i^ (6 ' 3) 

The integrand of (6.1) is finite at g = and takes value (2/T) ln(M/m) at this 
point. It is almost linear function and turns to zero at the upper limit q max . This 
means that radiative corrections ir rac i effectively screen the process at gluon momenta 
q < T (not gT\). The very integral can be roughly estimated at M >> m as 

~ 2 In — In — (6. 4) 

m IT y 1 

The result of numerical calculation is drawn at Fig. 4 by the bold solid line in com- 
parison with the thin solid line of the Born's term. In the numerical calculations 
of Ref. [8] no radiative corrections was introduced and the gluon momentum was cut 
off at q ~ gT. So the order of magnitude of the effect derived in this way may be 
different in the parametric scale (obviously there is no difference at g ~ 1). The finite 
part of the ir a bs is at least M/T times less than the leading term and the rest part 
of the radiative correction is suppressed by the extra factor exp(— M/T). They are 
plotted by the dashed lines. 

The induced emission along with the corresponding radiative correction has a 
statistical weight CgCq- The IR-finite combination of the leading terms from n em and 
7r ra d reads as 



<J - ^^f^ liGe-^M 2 J™ f e-^{-(l - ^)(L(q)e-^ - L(0)) + 

+ (l + 2 i£)(v(q)e^ T -v(0))} (6. 5) 

The integrand of (6.5) takes value (1/T) ln(M/m) at q = and the thermal 
exponent effectively cut off the upper limit at q ~ T. So the order of this integral is 
about — ln(M/m) which is also M/T times less than the leading term and produces 
negative contribution to the total rate of emission. One more term related to this 
process is positive and suppressed by the factor exp(— M/T). This curves are plotted 
at the Fig.5 by the thin solid lines. This contribution to the rate of emission is small. 
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The Compton process and the radiative corrections due to %^ ad do not give much also. 
They are plotted at Fig.5 by the dashed and dotted lines respectively, apart from the 
factors exp(— M/T) and Cq,g 

The spontaneous emission of a gluon by its statistical weight (q is going along 
with the vacuum vertex radiative correction n^ t . The IR-divergent part of the real 
process and that of vacuum vertex together give 



<r « ^^e-^M* f - ^)(me-^ - L(0)) + 



2m 2 

+ (1 + M2(v(q)e- 9/T ~ v(0))} (6. 6) 

It is easily estimated as — (2m /T) \n(M/m). Sign minus means the negative contribu- 
tion of this term to the rate of emission. The factor before log has replaced a standard 
ln(m/A) in course of IR-cancellation. This change is due to the shape of quark dis- 
tribution and we have no reasons to introduce an artificial mass of a photon. But the 
collinear singularity, the big negative — log 2 (M/m), originating from the nonsingular 
integral (4.15) in the Eq. (4.13) had survived. All these terms are plotted at Fig.6. 
The bold solid line presents (6.6). Thin dotted lines are due do the IR-finite terms of 
this process. 

The bold dashed line relates to the vacuum vertex. We notice that in Refs.[5,6] 
just this term has killed main contribution (6.4) from the process 3 — > 1. This is a 
remarkable consequence of the detail balance in the equilibrium heat bath. In our 
case this can not happen because these terms have different powers of fugacities. 

An appearance of the log 2 is a known phenomena and need the further summation 
of double logarithms. To understand what should be the result of this summation 
we remind that spontaneous emission of the gluons and the leading ladder diagrams 
of the vertex work together and result in the exponential K-factor. Related only to 
spontaneous processes, these diagrams have the common weight (q . We have already 
got a cancellation of the infrared singularities between soft real gluons and that in 
the vertex. So we may estimate the total yield of these two processes as Born's term 
times K-factor, 

a s C F 2 M 2 

fT~exp( — In — ), (6. 7) 

2n m z 

that is as strongly suppressed and even vanishing in the limit of m — 0. It is given 
by the thin dashed line at Fig.6 (with m ~ gT). If we had left the logarithmic term 
(4.15) as it is then we could get the negative total rate of emission in case of Cg < (q 

A spontaneous emission of soft gluons which naturally accompany the Coulomb 
interaction of the particles with the hard momenta are considered as the common 
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process with the very interaction of the particles. All parameters of the partons 
distribution are "external" with respect to this block and do not influence its internal 
dynamics in the gluon sector. It even does not seems logical to replace the mass m 
in the .fT-factor by the thermal-one. Though it came from the quark propagator with 
big momentum p° ~ |p| ~ M/2 (see Appendix B) the internal length scale of the 
processes convoluted to the i^-factor is 1/T which is less than the length scale 1/gT 
given by the forward scattering (see below). 

An exponentiating of the double-logs is in line with the accepted above scheme 
of renormalization of the vacuum parts of the self-energy and the vertex. Indeed, 
an assumption that the UV-divergence is removed by subtractions on the asymptotic 
mass shell means that we consider the interaction with the vacuum fluctuations as the 
dominant one and do not allow the interaction with the partons to destroy the mass 
shell. The existence of such a dominant part of interaction is a footing of any normal 
perturbation theory for the many-body systems. In fact, this is the main statement 
which allows to start with the initial distribution of the free partons. 



7. Conclusion and Discussion. 

This paper presents an analytic calculations of the rate of emission of heavy dilep- 
tons from nonequilibrium quark-gluon plasma adjusting the final result to the needs 
of computer simulation of the parton cascade. Calculations are based on the real-time 
kinetic approach developed earlier [1,7,11]. 

An explicit analytic calculation demonstrate that mass singularities remain in the 
final answer for the rate of the heavy dileptons production from the quark gluon 
plasma. We show that the previously obtained opposit result [5,6] is an artifact of 
the true thermal equilibrium in the partons background. One should not expect that 
in heavy ions collision this kind of equilibrium will be reached. This is in agreement 
with the known fact that at high momentum transfer the double logarithms in the 
physical answer is rather a rule than an exception [12]. 

The main practical consequence of this calculation is a conclusion that the main 
process beyond the first Born's one is annihilation of quarks with the absorption of a 
gluon. The kinematic of this process allows dilepton to store an additional amount of 
the internal energy. At M/T ~ 10 this process gives as much dileptons as the direct 
process. Ratio of its rate to the rate from the Born's term grows up almost linearly 
with the increasing of the dilepton mass, 

f/A ' ' ^ ;(n Born (k) + n a Uk)), (7. 1) 



d^kd^x 3(2tt) 6 M 2 

Aix 
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^ ie^N c C F 2 - W t M 2 1ti M 

The processes of annihilation with the spontaneous and induced emission of a gluon 
and the Compton process give a vanishing contribution in compare with this process. 

The a priori unexpected result of the above calculations is the cancellation of all 
those singularities in the thermal graphs which originate from soft gluons. This seems 
to be very significant as it reveals some hierarchy of scales and is worth of physical 
discussion. 

We got a new scale: the effective cancellation takes place up to gluon momenta 
~ T and this natural cut-off is greater than an obviously supposed gT. The later is not 
so surprising mathematically as we compare contributions of the graphs of the same 
perturbation order. Remember that the gluonic exponent in the finite combinations 
like (L(q)e~ q / T — L(0))/q have appeared in the Eqs.(6.1-5) only from the specific 
form of the quark distributions and from the conservation laws and without any 
connection with the distribution of the gluons in the partons cascade. Notice that 
if statistical weight in this combination was equal to 1 (which is the case of process 
in the free space) it would not contain the mass singularity also. This would have 
been in complete agreement with the KNL theorem which implies an absence af any 
statistical weights for the intermediate states. Probably the true thermal equilibrium 
is a unique case of populated phase space when a complete cancellation of the mass 
singularities takes place. 

The observed cancellation at q < T means that the two processes, spontaneous 
emission of a real gluon and an exchange by the space-like gluons, are mutually 
coherent at all gluon momenta q < T. As only a sum of these processes is physical it 
means that gluons with q <T do not participate this process at all : any attempt of 
the quark to emit the long-wave gluon is immediately killed by the rearrangement of 
the proper fields due to static interaction. This is in agreement with the typical size 
1/T 3 of the volume per parton which is free from the "third bodies". 

The collinear singularities did survive and they do need the additional cut-off. 
The later can be estimated from the physics of the parallel geometry of the front 
collision. In this geometry the gluon and the massless quark interact infinitely long 
until at least one interaction with the third parton will interrupt this process and 
destroy the mutual coherence of the real and virtual diagrams. 

The external cut-off should come from the smallest of the three lengths. The first 
one is the Compton wave-length l c ~ 1/m. The second one is the mean free path 
defined via the scattering cross-section, l m f p ~ \fg 2 T. The third one is connected 
with the amplitude of the forward scattering, lf s ~ 1/gT ~ l/m therm . As long as we 
keep quarks massless and consider coupling as the small parameter we must chose 
the length lf s as the smallest one. So whenever we meet the collinear singularity the 
momenta of both real and virtual quarks should be cut off from the below at ~ gT. 
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More formal arguments in favor of the scale gT come from calculation of the quark 
self-energy given in Appendix B. We could not use the previous results by Klimov and 
Weldon [13] as we deal with essentially different statistical background. It sufficiently 
changes (simplifies) an analytic structure of the quark self-energy. As it is clearly seen 
from equations (B.6), (B.7) and (B.ll) at high quark momentum the main deviation 
from the dispersion law of a free quark is formally simulated by the "thermal mass" 
of a quark, m t herm ~ gT. This does not mean that quark really acquire an additional 
rest mass. Indeed, the m 2 herm dominates in the dispersion equation (B.6) only when 
|p| >> T when its influence onto balance of the energy and momentum is negligible. 
We conclude that in the processes with hard quarks a thermal mass establishes the 
physical boundary of the possible coherence length in different processes. This is 
the main reason why we prefer to use a cumbersome "length associated with the 
amplitude of the forward scattering of the quark in partonic media" instead of a 
short and simple "thermal mass". 

An actual need of such a boundary appears only in collinear configuration of the 
hard processes when its geometric origin is quite evident. 

In the kinematic region of the heavy dileptons production (as well as for hard 
photons and low-mass dileptons of high energy) an approximation of the quark dis- 
persion law by formal introduction of thermal mass occurs to be remarkably exact. 
The next terms are suppressed by a small factor ~T/|p| ~ T/E. But this is not true 
for the case of the small time-like momentum of a quark when these next terms are 
very large and even divergent. 

A significant piece of work was done by R. Pisarski [14] in attempt to overcome 
specific problems of the soft region, |p| < gT, by means of selective resummation of 
the perturbation series. This extremely difficult study is confined to the ideal thermal 
equilibrium. Probably a final solution will demand a nonperturbative approach. 

For the heavy dileptons production when the dilepton mass exceeds any of other 
scales the particular choice of the cut-off scale is not so significant within the reason- 
able accuracy. The cut-off mass appears only under the logarithm and its variation 
is not so valuable. 

A choice of the cut-off may cause the qualitative changes in calculations of the 
emission of the low-mass dileptons when we are close to the threshold of the existence 
of the very Born's term along with the cascade part of the radiative corrections. They 
all appear only at M 2 > 4m 2 (~ 4g 2 T 2 ??) and are absent at the lower masses of the 
dilepton. 



I am indebted to G. Brown, E. Shuryak and the Nuclear Theory group at SUNY 
at Stony Brook for continuous support. 
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Appendix A. Singular integrals 

In course of calculations of terms associated with thermal parts of the quark self- 
energy we had met extremely singular integrals like 

/ = / d 4 pS(p 2 - m 2 )5((p - kf - m 2 )9(p )9(k - p ) J F(p) (A. 1) 
J p A — m A 



were F{p) is a regular function. We choose reference frame |k| = and get 

r I ,*? r- / 9 9\ r// , \ 9 9\ ' 

Po -w 

where oo 2 = p 2 + m 2 . Using the known formulae, 



= dp J d 3 p5(p 2 - uj 2 )5((p - k ) 2 - uj 2 )- -F( Po , p) (A. 2) 

JO J Vn — LU Z 



V I 1 11 

2- = — — + — — , 2m5{x) = —— + 



x x + ie x — ie x — ie x + ie 

Am5{x)- = - Krj + 7 j ry 2 (A. 3) 

x (x — ie) 1 (x + it) 1 

and absorbing all the regular part of the integrand into the function 

GM= ~, - f' P °' P) , r (A. 4) 

(p + uj) 2 {po-k -Lu) 

we get 

1= f d3p f G(p , p)dpo 

J 8tt 2 J(c) (p - uj) 2 (p -ho- to) 2 

contour (C) is given by (k — u) — ; u—; (k — Now it is easy to see that 

the internal contour integral differs from zero only at a single point u — k = k : 

I = - l -jd^G'{u^)5{u-^) (A. 6) 

This equation allows to perform the remaining integration over <i 3 p and being applied 
to Eqs.(4.3) and (4.4) results in Eqs.(4.10) and (4.11) for -K mass - 

Appendix B. Quark self-energy for a nonequilibrium partons 
distribution. 
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The invariant decomposition of the "thermal" part of the self-energy of the massive 
quark against the background of the initial partons distribution can be easily found, 

M r ^ = -^^{Am\+p[r ] -(pl + \p\ 2 )\-2po^]-u[poV-PoP 2 \-2 P ^]}, (B. 1) 
47r|p| z 

Where po = pu and |p| 2 = p\ — (pu) 2 . The extra |p|~ 2 in quark self-energy originates 
from the existence of the preferred reference frame and provides additional parameter 
like T/\p\ which is small in the region actual for heavy dileptons. In the rest frame 
of the media we will wright it in the form 

M£ t = W + 7°[/ + 7pV (B. 2) 



W = |^4m(A B + X F ), U = -2|^(£b + & + Po^b), (B. 3) 

a s C F riB + rip p 2 Q + |p| 2 - m 2 p 2 - |p| 2 + m 2 2p 
47T |p| 2 |p| 2 |p| 2 |p| 2 

with the following notations, 

2 f 2 r 

r) B = - d 4 q5(q 2 )n B (qo), rj F = - / d 4 qS(q 2 - m 2 )n F (q ), (B. 4) 

71 J 71 J 

Xb = - d q— n B (g ), A F = — / d g-^- ^ — ra F (g ), 

7T J p z — mr + 2gp 7r J + m z — Iqp 

£b = ~ d q— q nB{qo), £f = — / a q-r- — q n F {q ) 

7T J p A — m z + Iqp 7T J p z + m z — Iqp 

The dispersion equation, det [p — m — M ret ] = 0, can be rewritten in the next form, 

p 2 - |p| 2 - m 2 + 2(p U + \p\ 2 V - mW) + (U 2 - \p\ 2 V 2 - W 2 ) = 0. (B. 5) 

The last term of this equation is of the order a 2 . So if we are looking for the roots 
near their unperturbed positions, relaying upon the small value of a s , we have every 
reason to omit it. Up to the a s -order this equation is as 

p 2 - |p| 2 - m 2 - ^[rj B + Vf- (P 2 ~ Sm 2 )(\ B + X F )} = 0. (B. 6) 

In the case when the Lagrangian mass obey the inequality, m « « T « M, 
the explicit values of the invariants are as follows (E 1 is the integral exponent), 
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VB = 8( G T 2 , r) F = 8( Q T 2 (B. 7) 
\b = CgL( P ) \ f = CqL( P ), (B. 8) 

L(p) = l-[ e b- El ( b _) _ e ~b- E 1 (-b.) - e b +E 1 (b+) + e- b +E l {-b + )} 

IpI 

Z B = ( G K(p), £ F = ( Q K(p), (B. 9) 
Kip) = ^y-E^il - &_) + e- b - J E 1 (-6_)(l + 6_) - 

IpI 

-e^EMil - b + ) - e^E 1 (-6 + )(l + b + )} 

where b± = (po ± |p|)/2T and x) = E 1 (e~ l7T x) We are basically interested in the 
values of quark self-energy at M » T, p + ~ M and p_ ~ 0. In this region we can 
use an approximation, 

b + e b +E 1 (b + ) w 1 and ^(ft-) « -C B - ln(6_). (B. 10) 

It gives simple expressions, 

4T 2 2T T 2 o_ 

L (P)~-T^- OT T^' K{p)^—[2-2C E + m-2\n^-} (B. 11) 

IpI IpI IpI ^ 

So, as long as we are interested in corrections to the propagator specific for the heavy 
dileptons emission, when |p| ~ M/2, the term ~ a 2 in the Eq.(B.5) is suppressed 
by the additional powers of T/M. Up to the next orders in a scale of T/M only the 
invariants i] remain significant and we can replace Lagrangian mass by the thermal 
one, 

m 2 herm = ^(CQ + C G )T 2 (B. 12) 

Contrary to a simple calculation of the Green function for the field in QGP the 
problem of the field quantization in terms of quasiparticles may have no solution. 
Indeed, to be successive on this way we should expand the propagator using the 
complete set of the eigenf unctions of the wave equation modified by the self-energy 
and obeying a certain dispersion law. Even in classical plasma such a set does not 
exist. The Green function of the Landau kinetic equation allows to solve any problem 
with the given initial data but a complete set of the solutions of kinetic equation for 
longitudinal waves in a plasma form the so-called Van-Kampen modes which do not 
obey any dispersion law [15]. 
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